Unusual Transition Patterns in Bose-Einstein Condensation 
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We analyze the possible transition patterns exhibited by an efltective non-relativistic field model 
describing interacting binary homogeneous dilute Bose gases whose overall potential is repulsive. We 
evaluate the temperature dependence of all couplings and show that at intermediate temperatures 
the crossed interaction, which is allowed to be attractive, dominates, leading to smooth re-entrant 
phases. At higher temperatures this interaction suffers a sudden sign inversion leading to an abrupt 
discontinuous transition back to the normal gas phase. This situation may suggest an alternative way 
to observe collapsing and exploding condensates. Our results also suggest that such binary systems 
may offer the possibility of observing Bose-Einstein condensation at higher critical temperatures. 
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The study of symmetry breaking (SB) and symme- 
try restoration (SR) mechanisms have proven to be ex- 
tremely useful in the analysis of phenomena related to 
phase transitions in almost all branches of physics. Some 
topics of current interest which make extensive use of 
SB/SR mechanisms are topological defects formation in 
cosmology, the Higgs-Kibble mechanism in the standard 
model of elementary particles and the Bose-Einstein con- 
densation (BEC) in condensed matter physics. An al- 
most general rule that arises from those studies is that a 
symmetry which is broken at zero temperature should get 
restored as the temperature increases. Examples range 
from the traditional ferromagnet to the more up to date 
chiral symmetry breaking/restoration in QCD, with the 
transition pattern being the simplest one of going from 
the broken phase to the symmetric one. 

A counter-intuitive example may happen in multi-field 
scalar models, as first noticed by Weinberg iQ. Consid- 
ering an 0{N) X 0{N) invariant relativistic model, with 
two types of scalar fields and different types of self and 
crossed interactions, Weinberg has shown that it is possi- 
ble for the crossed coupling constant to be negative, while 
the model is still bounded from below, leading, for some 
parameter values, to an enhanced symmetry breaking ef- 
fect at high temperatures. This would predict that a 
symmetry which is broken at T = may not get restored 
at high temperatures, a phenomenon known as symmetry 
nonrestoration (SNR), or, in the opposite case, a symme- 
try that is unbroken at T = would become broken at 
high temperatures, a case called inverse symmetry break- 
ing (ISB). Since then the model has been re-investigated 
by many other authors using a variety of different meth- 
ods, both perturbative and nonperturbative analytical 
and numerical methods, giving further support to the 
idea of SNR/ISB. A recent review [|j lists most applica- 
tions and gives an introduction to the subject, discussing 
other contexts in which SNR/ISB can take place in con- 
nection with cosmology and condensed matter physics. 

Two of the present authors have also treated the prob- 
lem nonperturbatively taking full account of the cumber- 
some two- loop contributions M. The results obtained in 



Ref. ||3|] were shown to be in good agreement with those 
of Ref. [Q, where the SNR/ISB phenomena were studied 
using the Wilson Renormalization Group (WRG) and the 
explicit running of the (temperature) dependent coupling 
constants has been taken into account, showing that in 
fact the strength of all couplings increase with the tem- 
perature, enhancing SNR/ISB. This result completely 
rules out a possible decrease of the strength of the nega- 
tive crossed coupling that would lead to the eventual SR 
at higher temperatures. These interesting results from 
finite temperature quantum field theory raise important 
questions regarding a possible experimental observation 
of those phenomena in current laboratory experiments. 

The recent experimental achievement of BEC and its 
further improvements has opened the interesting possibil- 
ity of probing and studying finite temperature quantum 
field models and methods, currently used in cosmology 
and particle physics, in the laboratory. One of the re- 
markable things about BEC of dilute atomic gases is the 
possibility of adjusting several experimental features by 
fine-tuning the parameters with a high level of control 
and accuracy (for recent reviews see Ref. [|). We could, 
for instance, envisage the possibility of investigating the 
SNR/ISB phenomena using a system composed by a mix- 
ture of coupled atomic gases, like the ones recently pro- 
duced 1^ in which one has the same chemical element in 
two different hyperfine states and that may be treated as 
"effectively distinguishable" , or just consider the mixing 
of two different mono-atomic Bose gases. 

Here, we shall do a theoretical investigation of the 
possible transition patterns followed by this type of sys- 
tem paying special attention to the SNR issue since one 
expects, on basic grounds, that as the temperature in- 
creases the atoms will leave the condensate phase going 
to the symmetric phase of a normal gas. As a byproduct 
we shall see how the crossed couplings shift the critical 
temperatures with respect to those obtained for uncou- 
pled gases, creating the possibility of obtaining BEC at 
higher temperatures. 

The model we consider is similar to the ones used in 
other theoretical studies of homogeneous dilute coupled 
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Bose gases 0, that consists of a hard core sphere gas 
model described by non-relativistic interacting (complex) 
scalar fields, with an overall repulsive potential. This 
system can be described by the following U^{1) x U^{1) 
invariant finite temperature Euclidean spacetime action 



Se{I3)=J drjd^xi^* 
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where, in natural unities, T — 1/(3. The associated chem- 
ical potentials are represented by (i — tp or 0) while 
rrii represent the masses. For the hard core sphere self- 
interactions we take the phenomenological coupling con- 
stants as being the ones normally used in the absence 
of crossed interactions and which are valid in the di- 
lute gas approximation In terms of the correspond- 
ing s-wave scattering lengths they can be written as 
gi = Sirai/mi. To make contact with the analogous po- 
tential used in the prototype relativistic models for SNR, 
we take the overall potential as being repulsive, bounded 
from below. This requirement imposes the constraint 
condition > 0, > and g^g^ > -ig^. 

At T — 0, in the absence of crossed interactions 
and for /i^ > the classical potential of Eq. (|^) ex- 
actly reproduces the case of an uncoupled gas with sym- 
metry broken ground states representing "condensates" 
= {2fi^/g^y/^ and = {2iJ.^/ g^y/"^ . At the same 
time, the case fj,^ < and /i^ < corresponds to the nor- 
mal symmetric phases without any condensates, It/"! = 
and \(f>\ — 0. Our aim is to investigate, for non-negligible 
binary interactions, how the inclusion of temperature cor- 
rections can alter this picture. This can be achieved by 
computing the temperature dependent effective chemical 
potentials (see for instance Ref. |^ for an analogous anal- 
ysis for the mono-atomic gas) defined as a solution of the 
gap equation /Xi(T) = /Lti(O) -I- Ef (p), where Sf (p) is the 
field temperature dependent self-energy. 

The phase structure of the model is then given by 
the sign of HiiT) at a given temperature. One has 
/ii(T < T*) > in the broken condensate phase and 
/ii(T > Tl) < in the symmetric normal-gas phase. 
At the same time the Hugenholtz-Pines theorem imposes 
/ii(T = r*) = 0. Our next step is to evaluate the thermal 
self-energies which we shall do in a nonperturbative self- 
consistent fashion so as to avoid any potential problems 
associated with the usual perturbative calculation (see, 
e.g., Ref. |^ for a discussion in relativistic field theory). 
One can perform a one-loop self-consistent resummation 
by using the effective dressed propagator D,i^i.(p,w„) 
[—iujn + uJi]~^ where w„ = 2?™//? are the bosonic Mat- 
subara frequencies and uJi = p^/ (2mi) — ^ii{T). One gets 
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where gi{Q) and g{Q) indicate the bare, zero temperature, 
coupling constants. The sum in Eq. (||) can be easily 
performed and the resulting momentum integrals lead 
to well known Bose integrals For simplicity, in the 
following we consider atoms (fields) with approximately 
the same mass, ~ = m. One obtains two one- 
loop self-consistent coupled equations given by 



^i,{T) = ^,(0)-.g,(0) j Li3/2[exp(/3^,(T))] 

/ \ 3/2 

- 3(0) (^^j Li3/2[exp(/3Ai,(T))] , (3) 

where Li„(z) is the polylogarithmic function. As shown 
below, /ii(0) << Tc for realistic parameter values. We 
shall also see that ^i(T) quickly decreases (in the SR 
case) or remains approximately equal to fJ-i(0) (in the 
SNR case) with increasing T so that one may safely con- 
sider the high temperature approximation ^i{T)/T << 1 
by taking Li3/2[exp(/3/Xj (r))] ~ C(3/2). One then obtains 
the critical temperatures 
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Eq. (^ displays some unusual effects due to the presence 
of crossed interactions. There are three interesting cases 
which depend on the sign and magnitude of 17(0). Taking 
g{Q) > one observes a shift in the critical temperatures 
indicating that the BEG / normal-gas transition occurs at 
lower temperatures compared to the usual mono-atomic 
case (5(0) = 0). If 5(0) < but 1^(0)1 < g,{0) then the 
transition occurs at higher temperatures. Despite these 
important quantitative differences symmetry restoration 
does take place in both cases. Now consider the case 
where g{0) < but \g{0)\ > 5^(0) (in this case the bound- 
ness condition assures that |5(0)| < 5^(0)). For this situ- 
ation something unexpected occurs concerning the (j) field 
since Eq. (^ does not allow for a finite, positive real crit- 
ical temperature value. This is a manifestation of SNR 
within our two-field model and is analogous to what is 
seen in the relativistic case. At the same time the field ijj 
suffers the expected phase transition at a higher Tc com- 
pared to the g{0) = case. Obviously, which field will 
suffer SNR depends on our initial choice of couplings. 

For the following analysis let us set the parameters by 
choosing representative values for a mixture of gases such 
as ^^Rb and ^^Rb. Some realistic values are m = 86GeV 
— and a^f, = 2.5 x 10~'^(eV)^^, corresponding to the s- 
wave scattering length of ^^Rb, which fix the coupling 
g^{0) = 7.3(MeV)-2. Setting fj.^{0) = 5.67peV the value 
Tc ~ 32.5peV — 280nK is reproduced in the usual 5 = 
case. Throughout this Letter we set /^^(O) — fi^{0), 
^ keep 5^(0) fixed, while considering g^{0) and g{0) as tun- 

able parameters. In principle, this can be experimentally 
(2) achieved by appropriately setting magnetic fields close to 
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a Feshbach resonance as described in a recent application 
to ®^Rb As an illustration of the relevant phenom- 
ena we want to describe, we consider three possible sets 
of numerical values for (?0(O) and g(0), given, in unities of 
MeV"2 by: I) 5^(0) = 5.6 and 5(0) = 2.3; II) ^^(O) = 5.6 
and 5(0) = -2.3 and III) g^O) = 0.56 and ^0) = -0.84. 
Based on the general results given by Eqs. (|) and (|) we 
then have the usual transition patterns for both -0 and (f> 
for sets I and II where /i^ (T) and /i^ (T) go from positive 
to negative values, with transition temperatures that can 
be easily derived from Eq. (^). For set III, the field 
exhibits the usual transition from the condensed to the 
normal gas phase as the temperature is increased from 
zero. However, we have H^iT) > for arbitrarily large 
temperatures, which corresponds to SNR. This result is 
however misleading, as we demonstrate below. 

As a general result from the renormalization group, we 
expect that all coupling constants should run with the 
temperature For instance, if gi{T) decreases faster 
than g{T), or on the other way around, if the crossed 
coupling decreases faster than the self-couplings, the sim- 
ple analysis performed above can change drastically. As 
noted in the introduction, those substantial qualitative 
changes were not observed in the relativistic case. How- 
ever, a basic difference, in between the relativistic and 
non-relativistic model studied here, refers to the type of 
four-point functions allowed by each model. In fact, the 
contributions considered in the relativistic calculations 
include the t and u scattering channels as well as the 
s-channel contributions. At the same time, elastic and 
inelastic collisions are allowed whereas for the present 
case only elastic, s-channel contributions can be consid- 
ered since our effective model represents a system of hard 
core spheres. Keeping these facts in mind and perform- 
ing a nonperturbative one-loop calculation in terms of 
temperature dependent vertices and effective propaga- 
tors one obtains 



g(r)=.g(0) [l+g{0)B,Ak)]-\ 
where Bij (k) represents the bubble contribution. 
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The quantity Sij is a symmetry factor, s,;.j — 2 for i ^ j 
and Si J =5/2 for i = j, while flij is given by 



q2 + a(l - a)k2 
2m 



-a^^,iT)-{l-a)^lJ{T) , (7) 



where k represents the s-channel incoming momentum 
and a is a Feynman parameter introduced to merge the 
two internal bubble propagators in the derivation of Eq. 
(^. The q integral in Eq. is again performed with 
the help of the Bose integrals and one obtains 



Bjj(k) = s,j/3 
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da Lii/2[exp(-/3xij)], (8) 



where x^J = [a(l-a)k2/(2m) - /i,(T)a-^,(r)(l - a)] . 
Using the expansion Lii/2[exp(— z)] = 1.17 / z^/"^ — 1.46-1- 
0.2082Z - 0.0128z2 + 0{z^) (see Ref. §]) and setting 
k'^/2m = 3T (the average incoming two particle kinetic 
energy), the a integral in Eq. (H) can be performed. In 
the high temperature limit, one obtains 



S{T) ^ 9(0) 



3/2' 



(9) 



where — 3.7012. The equations for the couplings 
representing the self-interactions, gtp{T) and g4,{T), have 
exactly the same structure and can be obtained by replac- 
ing g with g^ or g^ in Eq. (||) and jij by 7^^^ = 4.6265. 
These results show that g^{T) and g,j){T) always de- 
crease, approaching zero monotonically at high temper- 
atures since g^{0) and (7</,(0) are positive quantities. The 
same happens for g{T), if g{0) is positive. However, 
for g{0) < 0, it can be easily checked that g^{T) and 
g^{T) decrease faster than g{T). Consequently the tran- 
sition picture described above will be shown to change 
in an unexpected way. In this case, Eq. (^) shows that 
g{T) is initially negative with an absolute which increases 
up to a certain critical inversion temperature given by 
Tg = (27r/m)3 [3.7012 5(0)]"^ where it develops a pole. 
Then, for T > Tg, g{T) suddenly becomes positive which 
means that, contrary to the relativistic case, one always 
has symmetry restoration at high temperatures in the 
nonrelativistic model of hard core spheres. It is impor- 
tant to note that this pole in Eq. (|^) has a completely 
different origin from the usual Landau pole found in the 
relativistic 0^ theory. Here the pole signals a phase tran- 
sition and if the nonperturbative self-consistent loop ex- 
pansion used here converges, as generically expected, a 
next order calculation will only lead to a small change in 
the value of Tg. 

It is instructive to look in more detail at sets II and III 
where ^(0) is negative. The results for HiiT) are shown 
in Fig. 1. Contrary to the previous analysis, with bare 
(temperature independent) coupling constants, we now 
have a completely unexpected behavior for both fields 
for some parameters values, as the ones exemplified by 
sets II and HI. For set HI, Fig. 1 shows that </> displays 
the typical SNR behavior while undergoes SR. The 
same figure shows that, for set II, both fields observe 
SR. In the three SR cases above, the transition from the 
BEC phase to the normal-gas phase is observed at critical 
temperatures between 380nK and 900nK. Then, as soon 
as the negative g{T) reaches a magnitude which is larger 
than that of gi(T), both type of atoms return to the 
BEC phase in a re-entrant transition. We note that, for 
set HI, the field ip also goes through a re-entrant phase 
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FIG. 1: The effective chemical potentials i-i^{T) (dotted and 
dot-dashed curves, for sets II and III, respectively) and fJ,4,{T) 
(continuous and dashed curves, for sets II and III, respec- 
tively) as a function of temperature (in units of eV). 
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FIG. 2: The high temperature discontinuous transitions suf- 
fered by both fields. Same notation as in Fig. 1. 

that happens at a higher temperature scale ^ 50/iK, not 
shown in Fig. 1. 

At the intermediate temperatures the boundness con- 
dition g^{T)g^{T) > 4g(r)2 does not hold, the over- 
all potential becomes highly attractive and unbounded 
from below. Then, suddenly, at the inversion tempera- 



ture T = Tg the potential becomes highly repulsive and 
both type of atoms leave the ground state at once. Con- 
trary to the first BEC/normal-gas transition, this second 
passage to the gas phase is a discontinuous transition due 
to the sudden sign inversion of g(T). This transition is 
shown in Fig. 2 for cases II and III for both ip and (p. For 
set II both fields go to the gas phase at a critical temper- 
ature Tg ~ 40/iK, while for set III, Tg - 340/xK. These 
values should be contrasted to the ones for the critical 
temperatures for the uncoupled gases, of T^ ~ 280nK, 
for both cases and T^ ~ 330nK and ~ ISOOnK, for cases 
II and III, respectively. 

Here, we cannot furnish more details concerning this 
discontinuous transition since, at the transition, our high 
temperature approximation and the effective model of 
two body interactions are not valid due to the very large 
values acquired by the fJ.i{T) and g(T). However, this 
should not invalidate the qualitative features of the tran- 
sition, at least on its neighborhood. 

In conclusion our results show that a system of cou- 
pled Bose gases with an overall repulsive potential, 
but attractive crossed interactions, may display unusual 
transition patterns, when the important temperature 
dependence of the couplings are taken into account. 
These patterns include the usual continuous transition 
condensate/normal-gas followed, at intermediate temper- 
atures, by an unexpected re-entrant, continuous, normal- 
gas/condensate transition. Higher temperatures induce 
a sudden change in the sign of the crossed coupling fol- 
lowed by a dramatic discontinuous transition of the type 
condensate/normal-gas. Those phases suggest a possible 
pattern for the observation of collapsing and exploding 
condensates. It is important to note that experimen- 
tally one knows that a change on the sign of the cou- 
plings may be achieved by adjusting external magnetic 
fields [0. However, our results further indicate that, 
at least for coupled systems, the temperature can also 
act as the external agent which drives the sign inversion. 
In addition, our results explicitly demonstrate that, con- 
trary to the relativistic case, and according to intuition, 
SNR cannot happen in the non-relativistic model of hard 
core spheres with temperature dependent couplings. Fi- 
nally, the coupled Bose-Einstein gas system with attrac- 
tive crossed interaction also seems to offer an avenue to 
obtain condensation at higher temperatures. 
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